Nonlinear a-Effect in Dynamo Theory 

by 

George B. Field, Eric G. Blackman, and Hongsong Chou 

Received ; accepted 



- 2 - 
ABSTRACT 

The standard two-scale theory of the dynamo coefficient a in incompressible 
isotropic helical MHD turbulence is extended to include nonlinear effects of B, 
the large-scale magnetic field. We express a in terms of statistical quantities 
that can be calculated from numerical simulations of the case B = 0. For large 
magnetic Reynolds numbers our formula agrees approximately with that of 
Kraichnan (1979), but disagrees with that of Cattaneo and Hughes (1996). 
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1. Introduction 

Magnetic fields of galaxies are important in astrophysics and cosmology. In astrophysics, 
because they enable fast particles to be accelerated and trapped, and affect the dynamics 
of star formation. In cosmology, because if galactic magnetic fields do not originate in the 
modern era, they could be relics from the early universe, carrying information about that 
period. 

Since Parker's (1955) paper on the a-fl dynamo, and his application of dynamo theory 
to the Galaxy (1971), most workers have attributed the origin of the magnetic fields of disk 
galaxies to the operation of an a-Q turbulent dynamo (Ruzmaikin, Shukurov & Sokoloff 
1988). However, the standard theory of the dynamo has always been open to criticism. 
For example, Piddington (1970, 1972abc, 1975ab) argued that the small-scale magnetic 
field produced by the small-scale turbulence required by the theory would rapidly grow to 
equipartition, quenching dynamo action, and this point has been demonstrated numerically 
by Kulsrud and Anderson (1992). 

Recently a series of papers (Vainshtein & Rosner 1991; Cattaneo & Vainshtein 1991; 
Vainshtein and Cattaneo 1992; Tao, Cattaneo & Vainshtein 1993; Cattaneo 1994; and 
Cattaneo & Hughes 1996; see also Seehafer 1994, 1995) have argued that dynamo action 
is quenched if the large-scale magnetic field B in velocity units exceeds a critical value, 

1/2 

B c = R M ' v , where Rm = v L/X is the magnetic Reynolds number of the turbulence, 
v is the turbulent velocity at the outer scale of the turbulence L, and A = i]c 2 /4:Tt is 
the magnetic diffusivity, with r\ the resistivity. They argue that this result is supported 
by direct numerical simulations of MHD incompressible turbulence (Cattaneo & Hughes 
1996). If they are correct, the classical a-Vt dynamo theory based upon Parker's (1955) 
paper and developed by the Potsdam group (Krause and Radler 1980; see also Moffatt 1978 
and Parker 1979), which applies to weak large-scale magnetic fields, is not applicable to 
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present-day galaxies, as B is observed to greatly exceed R M vq- 

The argument by Cattaneo and his collaborators depends upon their analysis of 
nonlinear interactions which are surely present. The effects of such interactions have also 
been considered independently in a series of papers on MHD turbulence by Pouquet and 
her collaborators (Frisch, Pouquet, Leorat, & Mazure 1975; Pouquet, Frisch, & Leorat 1976; 
Pouquet & Patterson 1978; Meneguzzi, Frisch & Pouquet 1981) using spectral methods. By 
and large this work supports the applicability of the classical theory even for large-scale 
fields B approaching t> , which encompass those observed. Chandran (1996) confirms the 
results of Pouquet, Frisch, & Leorat (1976) using a different spectral method. 

Since Rm ^ 1 in most astrophysical situations (as large as 10 20 in the interstellar 
medium of our Galaxy) the difference between the critical field advocated by Cattaneo 

— —1/2 

and his collaborators, B c ~ R M t> , and that implied by Pouquet and her collaborators, 
B c ~ v , is crucial. 

In this paper we extend the calculation of one of the dynamo coefficients, a, of the 
classical theory to include arbitrary values of B. Our result contains two terms, one which 
resembles that advocated by Cattaneo and his collaborators in its dependence upon Rm- 
The other resembles a formula proposed by Kraichnan (1979) on the basis of a simple model 
incorporating damping by nonlinear interactions. For R M ^> 1, the latter term dominates, 
so we find that the classical result for a applies for B of the order of v . We plan to extend 
the simulations of Cattaneo and his collaborators to other values of Rm in order to clarify 
our disagreement with them. 

We confine discussion to the alpha effect, and omit any discussion of turbulent diffusion, 
which is also controversial. 

Before proceeding to the nonlinear theory we review the classical linear theory. The 
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classical theory is based on a clear separation of scales between the scale L of the dominant 
turbulent motions (here called the "outer scale"), and the size D of the system, with 
L <C D. The large-scale magnetic field B satisfies an induction equation, 



where V is the large-scale velocity field (differential rotation in the case of a galaxy), and 
(v x b), called the "turbulent emf," is the spatial average of the cross product of the 
small-scale velocity v and the small-scale magnetic field b over a scale much smaller than 
D, but much larger than L. 

The first term in (1) leads to the so-called "Q-effect," according to which lines of 
force of B are stretched by the differential rotation, creating a growing toroidal field from 
a poloidal one. This term is not controversial, and will not be discussed further here. 
Controversy centers on the turbulent emf, (v x b) , which enables a growing toroidal field to 
feed back into the poloidal direction, giving exponential amplification of the large-scale B. 

The evaluation of (v x b) is usually restricted to incompressible isotropic turbulence. 
On might assume that if the turbulent velocity field v is isotropic, the small-scale magnetic 
field b would be isotropic also, so that the c-component of the emf would be 



in the light of an identity for isotropic tensors (Krause & Radler 1980). However, this is 
not correct, because B, being anisotropic, induces a term in b that is not isotropic even 
though v is isotropic. Evaluating this term requires the induction equation for b, which 
can be found by writing the induction equation for the total magnetic field B = B + b and 
separating off the small-scale parts to give 



dtb = V x (v x B) + V x (V x b) + V x (v x b) - V x (v x b) + AV 2 b . (3) 



t B = V x (V x B) + V x (v x b) + AV 2 b , 



(1) 



(v x b) c = e cde (v d b e ) = \e cde 5 de (v • b) = , 



(2) 
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The second term on the right in (3) represents a change of the reference frame to that 
moving at the mean velocity V, in which v may be assumed to be isotropic; we assume 
V = const., so we make that change and henceforth omit the term. The fourth term on the 
right is a large-scale quantity, so its time integral, being large scale, will contribute nothing 
when crossed wtih v (small-scale) and averaged, and so is neglected henceforth. The first 
term is — v ■ VB + B ■ Vv. As the v • VB term ultimately leads to turbulent diffusion of B, 
which is not our main interest here, we omit it, leaving 

d t h = B • Vv + V x (v x b) + AV 2 b . (4) 

The second term in (4) is neglected in the classical discussions in a step referred to as 
the first order smoothing approximation, or FOSA. Krause & Radler (1980) showed that 
this is legitimate if A is large, so Rm is small, because then Ohmic diffusion keeps b small, 
and the first term dominates the second. However, because Rm is large in astrophysics, this 
case is not relevant here. They also discuss the case when the Strouhal number 

s = x < 5 > 

is small, where r is the correlation time for eddies at the outer scale of the turbulence. 
Because L/vq = t e d, the eddy turnover time, S = r/t e d, and one might suppose that 
S = 0(1). In fact, it is observed experimentally that S ~ 0.2 — 0.3 in ordinary hydrodynamic 
turbulence (Pope 1994). Although this value is not as small as one would like, it provides 
a way to approximate (4). To understand how, it is important to distinguish between b^, 
the small-scale magnetic field when B = 0, and the perturbation to b when a small 
B is present. If the turbulent velocity field v is isotropic, b^ ^ is, and as shown in (2), its 
contribution to (v x b) vanishes. However, b^ is not isotropic, and its contribution to 
(v x b) does not vanish. 
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If we write b = b^ ^ + b^, is governed by 

d t b^ = V x (vW x b<°>) + AVV°) , (6) 

where v^ ^ is the isotropic turbulent velocity. Parker (1979; p. 511) shows that if A is 
small (the astrophysical case), @ leads to exponential growth of with a time constant 
t/S 2 until (in velocity units) it begins to approach i>' '. At that point, we expect that 
back reaction due to the Lorentz force associated wtih the large value of the small-scale 
field will result in a steady state in which the energy driving the turbulence at the 
outer scale (buoyancy forces in stars, supernova explosions in galaxies) is balanced by a 
nonlinear cascade to smaller scales, where it is dissipated by viscosity and/or Joule heating. 
This expectation is confirmed by the Pouquet et al. (1976) and by direct simulations, e.g. 
Cattaneo and Vainshtein (1991), Cattaneo (1994), and Cattaneo and Hughes (1996). As 
we explain later, we will take this steady state as the base state which is perturbed by B. 
The first term in (|G|) mediates exchange of magnetic energy with kinetic energy in MHD 
turbulence, as discussed in Appendix A, and it certainly cannot be neglected. 

However, we now argue that the corresponding term in the equation for can be 
neglected, as follows. Evidently is governed by the parts of (4) which are of first order 
in B, namely, 

dtb® = B ■ Vv<°» + V x (v<°> x b«) + AV 2 b« . (7) 

In the rest of this section, we assume that Rm 3> 1, so that we may neglect the third term. 
We formally integrate (0) to get 

b«(x,t) = B-vT d*iv^(x,ti) 

J — oo 

+ f dt x V x [v (0) (x,t!) x b«(x,ii)] . (8) 

This does not seem useful, because the desired quantity, b^, appears under the integral. 
However, we can show that the integral is much smaller than b^ itself, which appears on 
the left-hand side, if S -C 1. 
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Since is a stochastic function of t, only those parts of the integral in @ which come 
from times t\ which differ from t by less than a correlation time r will correlate significantly 
with v^ ) in the turbulent emf, /v^(x, t) x b^(x, t)Y Hence we can replace the lower limit 
on the integral by t — r, and then estimate the integral by r times the integrand at t\ = t. 
Since the order of magnitude of V is L -1 , the order of magnitude of the second term in @ 
is 

^bV{t) = SbV{t), (9) 

whose ratio to the magnitude of the left hand side of (||) is S. If, as we shall argue later, 
S < 1, we can neglect the second term in (||), and hence in (0), so, still neglecting diffusion, 
(0) becomes 

<%b« = 5 • Vv (0) . (10) 

Effectively ~ SB, so if 3 < 1, < B, even when B v(°\ However this says nothing 
about b(°\ which as we have stressed, approaches v ^' in value, even if > B. Hopefully 
this discussion clarifies a point which has led to confusion in the past. 

The solution to (|H]) can be written in component form as 

b^\t) = B p d p f Mvf\t'), (11) 

J — oo 

so that 

<V°) x b«) c = e cde (vf{t)B p d p fjt'vf\t')^ . (12) 

Averaging commutes with integration and differentiation, spatial differentiation commutes 
with time integration, and B is independent of time on the scale r, so 

(V x b) c = B p e cde f dt' (vf\t)d p vi°\t')) . (13) 

Since v^ ^ correlates with its derivatives for a time of order r, the integral is of order r 
times the average taken at a given time (Krause & Radler 1980). To make further progress 
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in what follows, we need an explicit dependence of the indicated correlation on t and t'. 
Because the turbulence is assumed to be steady, it can depend only on t — t', so a convenient 
representation is 

(v®mv^)=(v?>a,v? s >)e*^VT i (14 ) 

where the common time argument in the second average has been omitted because the 
turbulence is presumed steady. Hence (v x b) is given by 

(vxb) c = 5 p r ecde (ftf) . (15) 

At this point we use the assumption that the velocity v^ ) is distributed isotropically. 
According to Krause & Radler (1980), a third-rank isotropic tensor like that in fllBT) can be 
written 

^^^^(v^Vxv^)) , (16) 

so from (|13|), 

(v x b) = oB , (17) 

where 

a = _| r ^ v (o). Vx v (o)^ (18) 

is the classical expression for the dynamo coefficient (Krause & Radler 1980, Moffatt 1978), 
but with the additional feature that the turbulent velocities indicated refer to the zero-order 
state. The quantity in angular brackets, a pseudoscalar, is known as the kinetic helicity of 
the turbulent flow v(°). 

In an independent development, Pouquet, Frisch & Leorat (1976) calculated magnetic 
energy spectra for MHD turbulence, solving the spectral equations using a closure method 
known as the EDQNM (Eddy-Damped Quasi-Normal Markovian) approximation. They 
found that if (v^ ■ V x v^ ^ vanishes, the magnetic energy spectrum E^ 1 peaks near 
ko, the wave number at which turbulent energy is injected, and reaches a steady state in 
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which the total energy E^, the sum of the kinetic energy El and the magnetic energy 
Ejf 7 cascades to higher wave numbers, ultimately to be lost to Ohmic dissipation and/or 
viscosity. 

If, on the other hand, (v^ ■ V x v^y 7^ 0, they found that E^f inverse cascades, 
accumulating at an ever- decreasing wave number. They attributed this to the turbulent 
dynamo effect described above, operating in the nonlinear regime. In the special case that 
the wave number k of interest is <C ko, so that there is a clear separation of scales, they find 
an approximate expression for a quantity a k governing the growth of Ejf , where 

«*=-!/ dqe kqq (H^-q 2 H^). (19) 

J k/a 

Here a is a small parameter, 9k qq is an effective correlation time for modes of wave number q, 
and HY and q 2 E q M are the spectra corresponding to the kinetic helicity correlation function 
^v^(x) ■ V x v(°)(x + £)^ and the current (Keinigs 1983) helicity correlation function 
^b^(x) • V x b^(x + £)^>, respectively. (Here b is in velocity units, obtained by dividing 
b by \/4:iip, so that b is the vector Alfven velocity.) The term adopted by Pouquet et al. 
for (|19|), the residual torsality, has not reappeared in the literature. Here we note that the 
first term in (|i"9[ ) is similar to (|18D , so that classical dynamo theory can be interpreted as 
an inverse cascade in helical turbulence. The second term does not appear in the classical 
result, but Montgomery and Chen (1984) verified it using calculations in real rather than 
k space. Pouquet et al. thus found that dynamo action takes place in the fully nonlinear 
regime, with no restriction as to the magnitude of b or B, and that the classical expression 
for a must be modified by the addition of the second term in (|I9"D . 



2. The Method of This Paper 



In this paper we extend the classical analysis of the a effect into the regime of large B. 
We use the standard two-scale approximation in real space, and at several points we refer 
to the results obtained by Pouquet et al. (1976) using the spectral closure method. We find 
that when a clear distinction is made between the quantities v and b on the one hand, and 
their values when B = 0, and b^°\ nonlinearity due to B appears in a straightforward 
manner. A key feature of our derivation is the assumption that correlations are damped by 
nonlinear interactions, following Pouquet et al. (1976) and Kraichnan (1979). As shown by 
Kichanitov (1985), such an assumption can be justified by application of renormalization 
group methods. 

Our results depend upon the damping rate at each wave number k, 7^, the spatial 
spectra of v<W and b^, the value of B, and the value of A, the magnetic diffusivity, 
expressed in terms of the magnetic Reynolds number R M , which for convenience we assume 
to be equal to the Reynolds number R. We find that a is a well-behaved function of the 
parameters, and that the behavior for larger Rm is similar to that predicted by Kraichnan 
(1979). 

We take as our base state fully-developed MHD turbulence driven by external forces 
and in a steady state as a result of a turbulent cascade to large wave numbers, but with 
B = 0. As demonstrated by Pouquet et al. (1976), in such turbulence there is approximate 
equipartition between the magnetic energy Ejf and the kinetic energy E\ for wave numbers 
k~Z 3k , where L = k^ 1 is the outer scale of the turbulence. As shown by Parker (1979, 
p. 513), approach to the steady state occurs on the scale of the eddy turnover time 
t ed = L/v(°\ which, as stressed by Kulsrud & Anderson (1992), is much shorter than the 
dynamo growth time. According to the calculations of Pouquet et al, the saturation which 
occurs at small scales (large k) does not prevent the increase of magnetic energy on scales 
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larger than L if the turbulence is helical. It is important to note that the growth of a 
large-scale field as a consequence of an a effect does not substantially modify the spectra of 
Eff calculated by Pouquet et al. (1976) for k of the order of ko (their Fig. 8), so our concept 
of a base state independent of B is valid. We assume that the properties of the base state 
can be calculated once and for all by numerical simulations. Our results then allow us to 
calculate a in terms of those properties. 

Following Montgomery and Chen (1984), we present our calculations in terms of the 
Elssaser variables z ± = \ ±h, (Biskamp 1993) and have checked our results by carrying 
out the calculation in terms of v and b. Elssaser variables are naturally adapted to the 
problem, shortening the calculation substantially. More important, we show in Appendix A 
that unlike the kinetic energy E\ associated with v and the magnetic energy Ejf associated 
with b, both E£ and cascade directly in isotropic turbulence, allowing us to employ a 
single damping constant Ik = Ik = Ik t° describe the effects of the nonlinear terms in the 
base state. 

Our goal is to calculate the turbulent emf, (v x b). Since 



v = o 



i(z + + z-) 



(20) 



and 



b=|(z + -z-), 



(21) 



we have 



(v x b) c 



I((z+ + z-)x(z+-z-)) 



X z 



) 



c 




(22) 



We use this formula in what follows. 
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3. Evolution of the Turbulent Fields 

As explained above, the correlation indicated in (^) vanishes for isotropic turbulence, 
but is nonzero when one takes into account the perturbations of z + and z~ which are caused 
by B. To obtain these, we consider the dynamical equations for z ± , or equivalently, v and 
b. 

When Ohmic dissipation is included, the induction equation for b is, from (3), 

d t b = -v ■ VB + B ■ Vv - V • Vb + b ■ VV - v ■ Vb + b ■ Vv + AV 2 b , (23) 

where A is the magnetic diffusivity. As explained in §1, adopting in a frame of reference 
moving with V eliminates the third term on the right, and the fourth term can be neglected 
with respect to the sixth because we assume that the size of the system S ^> L. The 
remaining terms can be written in terms of B = B + b as 

d t b = -v • VB + B • Vv + AV 2 b . (24) 

The classical theory ignores the effect of B on the velocity, on the grounds that the 
Lorentz force associated with B is of order B 2 , hence negligible in the limit B — ► 0. 
However, as we have explained above, grows quickly to approximate equipartition, so 
even in a first-order calculation, a Lorentz force proportional to Bb^ must be included. It 
is therefore essential to consider the effect of B on v; we shall do so to all orders in B. To 
do this, we use the equation of motion for the small-scale velocity v: 

dtv = -v • Vv - Vp + B • VB - \7\B 2 + iA7 2 v + f , (25) 

where f is the applied force per unit mass and v is the kinematic viscosity; B, the magnetic 
field ~^Airp, is in velocity units. Hence 

d t v — -v • Vv - VP + B • VB + zA7 2 v + f , (26) 



where 



2 

We define the Elssaser variables for the field B as 
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P = p + \B 2 . (27) 



so that 



and 



Z ± = v±B = v±b±B = z ± ±B, (28) 



v = ±(Z+ + Z-) (29) 



B = b + B = i(Z+ - Z~) . (30) 



Then fl24T) becomes 



|9 t Z+-|5 t Z- = -I(Z + + Z-)- V(Z+-Z~) 

+i(z + -z-) • v(z+ + z~) 

+|AV 2 (Z+-Z-) . (31) 

(Note that we have set 9jB = because B varies only on the long time scale.) Equation 
(f26|) becomes 

|5 t Z+ + |5 t Z- = -i(Z+ + Z-)-V(Z+ + Z-) 

+i(Z + -Z-)-V(Z+-Z-) 

+ |z/V 2 (Z + + Z-) - VP + f . (32) 



Adding and subtracting (|31|) and ( p2|) yields 



9 t Z ± = -Z T ■ VZ± + + A)V 2 Z± + - A)V 2 Z T - VP + f . (33) 

The dissipative coupling between Z + and Z~ vanishes in the special case that the magnetic 
Prandtl number \jv equals unity For simplicity we assume that is the case in what follows; 
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we note that both Pouquet et al. (1976) and Cattaneo and Hughes (1996) also made the 
same assumption. Chou and Fish (1998) discuss the case Xjv ^ 1. 

From ( pffi ) the nonlinear term in ([33|) is 

Z T -VZ ± = (z T =p B) ■ V(z ± ± B) 

= z T • Vz ± t B • Vz ± , (34) 

where we have neglected VB for the reasons given previously. Hence (|33|) becomes 

d t 7^ = -z T ■ Vz ± + AVV - VP + f ± B • Vz ± , (35) 

where we have neglected d t ~B for the reason given previously. 

We adopt a perturbation expansion in B, in which the zero-order variables z^ - 1 
describe the turbulence exactly if B is zero. Therefore z ±( - ^ satisfies 

d t z^ = -z*°> ■ Vz ± (°) + AV 2 z ± (°) - VP (0) + f , (36) 

where 

p(0) = p (0) + | S 2|_ =o = p(0 ) + I (6 (0) )2 (3?) 

If we assume that the driving force is independent of B, f has the same value in both the 
zero-order and the perturbed state. 

Equation ( |36"D can be solved numerically for a variety of initial conditions, and so in 
principle any averages of zero-order quantities required can be computed. Our task is then 
to compute (v x b), which depends on the perturbations induced by B, in terms of averages 
over zero-order quantities. 

We let 

z± = z ± (°) + z'± , (38) 
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where z f± contains perturbations of all orders in B. Although in principle z f± can be 
represented as a power series in B, we find that it is not necessary to do so for the special 
case A/V = 1 . If we eliminate f between fl35[) and (|36|) , we find that 

(d t - XV 2 )z r± = -z'T . Vz ± (°) - z=F(°) ■ Vz'± - z^ ■ Vz 7± 

-V(P-P (0) )±B- Vz ± . (39) 



Inspecting (|39|) , we see that the left hand side and the term B • Vz f± on the right 



are linear in z /=t , hence easy to deal with. If Rm is moderately large, we can ignore the 
term AV 2 z' ± in assessing the order of magnitude of the remaining terms. Putting aside the 
pressure term for the moment, we can write (|3"9"|) in the symbolic form 



d t z' 



z'z^Az')\Bz^\Bz' 



(40) 



where the commas separate terms of potentially different orders, and where we have 
combined the first two terms on the right of ([39]) into one. As in our discussion of the 
classical case, we take Az' = J dtd z z', the change in z' after one correlation time r, to equal 
z' in order of magnitude, by definition. Then fl40|) implies that 



T 

ted 



z',(z') 2 /z (0 \B,Bz'/z^} , (41) 



where t ed = L/z^°> ~ Lj\ 



In the classical discussion it is assumed that because B is small, back reaction by the 
Lorentz force can be neglected, so the motions are hydrodynamic in nature. That allowed 
us to use the fact that r/t e d is small (see below). That is not really true, because, as we 
discussed earlier, the zero-order state quickly approaches equipartition. We want to argue, 
however, that even in the nonlinear case r/t e d is a small number. To do so, we appeal to 
the nonlinear calculations of Pouquet et al. (1976) regarding the zero-order state. They 
showed that the magnetic energy spectrum E^ 1 is in equipartition with the kinetic energy 



-17- 



spectrum E\ for 3ko, where fco is the wave number at which kinetic energy is being 
injected (the outer scale of the turbulence). In this range, back reaction is a major effect, 
causing the turbulence to become a field of interacting Alfven waves. However, these large 
wave numbers need not concern us, because in such waves, the current is perpendicular to 
the small-scale field, and the vorticity is perpendicular to the small-scale velocity, so that 
the two helicities vanish. 

Of crucial importance for us, Pouquet et al. (1976) found that in the range ^~ 3A;o, 
Ejf < E%, so that the motions are largely hydrodynamic in character, with only a modest 
back reaction of the magnetic field. This range is important because it contains most of 
the energy. Pope (1994) states that in a pure hydrodynamic turbulence it is experimentally 
observed that 

— = 0.2-0.3. (42) 

ted 

We assume that ( f£2|) applies to the energy-containing eddies in MHD turbulence driven like 
that of Pouquet et al. Although r/t e d is not a very small number, we may take it to be a 
small parameter e in analyzing (f|l]). Thus, ( f4~T| ) becomes 



z' ~ e 



z',{z') 2 /z^,B,Bz'/z^] . (43) 



It seems reasonable to assume that the balance in (]4"3T) is between z' and the third term on 
the right, so 

z ~ efi . (44) 

We check this by evaluating the first two terms under that assumption. We see that the 
ratio of the first to the third terms, and the ratio of the second to the fourth term, are both 
e, so the first and second terms can be neglected as a first approximation. Note, however, 
that the ratio of the fourth to the third term is 

z'/z {0) = eB/z {0) ~ eB/v {0) , (45) 
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which cannot be neglected because, although e is small, we want a result which is valid to 
all orders in B/v^°\ Indeed, this is the source of the nonlinearity in our calculation. 

Note that although we have used the short correlation time approximation to simplify 
our equations for z', no assumption is made regarding the ratio of b to B. This is a 
step forward, because the classical discussion has been justly criticized for ignoring back 
reaction, which amounts to assuming that b is small. 

We now apply the divergence operator V- to (|39"D with the first three terms on the 
right dropped. Because B is constant and V- commutes with B • V, the fact that 

V • z ± = V ■ (v ± b) = (46) 

for incompressible turbulence then implies that 

V 2 (P - P {0) ) = , (47) 

which for a homogeneous system implies that P — P^ = const., so that 

V(P - P (0) ) = . (48) 

Hence (p9| ) becomes 

(d t — AV 2 =F B • V)z' ± = ±B ■ Vz ±(0) . (49) 
We must solve this equation in order to calculate the turbulent emf according to 

The operator B • V is best handled by introducing the spatial Fourier transform 

z'^k, t) = (2tt)- 3 J dxe- ik - x z /± (x ) t) , (50) 

so that when written in terms of the components of z, fl49|) becomes 

(d t + AA; 2 ^2k-B)^ ± (k,t) = ±tk-Bzj i0) (k,t) , (51) 
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to which the solution is 



zf{k,t) = ±zk-B / dtiexp{- [^ik-B + Xk 2 ] (t - tj)} zp 0) (k, £ x ) 



(52) 



where we have put the lower limit equal to — oo because we will find that the short 
correlation time makes values of t\ significantly smaller than t irrelevant. The nonlinear 
dependence on B is evident here. 

From ( ^2] ) and (|3"8"D, the term of zero order in z^ drops out according to (2), leaving 

+(0) 



vxb > c = -\(e cd ez7 \x,t)z'-(x,t))-U6 cde z' d + (x,t)z;(°\x,t) 



-h (e cde z' d + (x,t)z' e (x,t) 



(53) 



Inverting the Fourier transforms, we see from ( p0|) that 



±(o) 
z d 



;x,t) = j dk^z^\k,t) 



(54) 



^ik'-B + AA;' 2 ] (t-t 1 )}^ (0) (k / ,t 1 ) . (55) 



and, from fl52|) 

z^(x, t)=±ij dk'e ik ' x k' ■ B j dh exp { 
Hence the first term in (|53| ) is 

liecde/J dk J rfkV (k+k ' )x k • B |* tZtiexp{-[ik'.B + AA;' 2 ](t-ti)}4 (0) (k,t)2-^(k',t 1 

(56) 

Because averaging commutes with integration, this can be written 

lijdkj rfkV (k+k,) - x k'-B y* dhexp{-\ik! • B + AA;' 2 ](t -t x )} (e C(ie z+ (0) (k, t)z e ~ (0) (k', *i 



(57) 



The indicated correlation decreases rapidly to zero as k — k and t — ti go to zero. In 
Appendix B we show that the following representation is a reasonable one: 



e„ te 4 (0) (k,t)l e -W(k',*i)) = e cde ^ e (-k)5(k + k / )e-^l*-' 1 l . 



(58) 
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In principle, 7^ could be different for z + and z~ modes, but we show in Appendix A that 
in turbulence in which the cross helicity K — | (v • B), an ideal invariant, vanishes, the 
associated energies E + and E~ cascade directly and are equal. It is therefore reasonable 
to assume that if K = 0, 7^ = 7^. Following Pouquet et al. we assume that K — 0, and 
thus, that 7^ = 7/7. Note that because 7^ is a parameter of the zero-order state, it does not 
depend upon B. From Krause and Radler (1980, p. 75) an isotropic tensor like R de (k) can 
be expressed in the form 

R de (k) = (A-k- 1 d k B)5 de -(d 2 k B-k- 1 d k B)k- 2 k d k e 

+ ik~ x d k Ct deS k s , (59) 

where A, B, and C are functions of k alone. Therefore 

e cde R de {-k) = -2ik c k- l d k C (60) 

is nonvanishing only for turbulence lacking reflection symmetry, and so ( p^ ) becomes 

;e cde 4 (0) (M)C (0) (k',ti)) = -UKk-^Ce-^^ . (61) 



The integration of (|57|) over k' converts k' to — k, so (k + k') - x = 0. Since \t — ii| = t — t\ 
in the integration over t±, we find that fl5"T|) equals 

- j dkk- Bkck^dkC J dt x exp |[zk • B — \k 2 - 7fc ] (t - t t ) } 



We note that the second term in ( |5Tf ) can be written 

+ ±e cde (<+(x,t)z7 (0) (x,t)\ (63) 



because of the e cde . Then reference to (|54]) and (|55D shows that this is the negative of (]6 
with i <-» —7 Hence the sum of the first two terms in (|53|) can be written 

ik • B ik B 



z y dkk c k 1 d k C 



(\k 2 + 7fc ) - ik ■ B (A/c 2 + 7fc ) + zk ■ B 
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-2 y dkkck^dkCk ■ B(Xk 2 + -f k )[{Xk 2 + 7fc ) 2 + (k ■ B) 2 ]" 1 . (64) 



We now turn to the third term in (p3D, which is 

-^JdkJ rfk'e'( k+k ' ) x (k-B)(k'-B) 
■f diiexp{-[-zk- B + Xk 2 ](t - tM 

J — oo ' 

■ f dhexp {-[ik' ■ B + A£; /2 ](t - t 2 )} (e crfe l+ (0) (k,t 1 )4-(°)(k , ,t 2 ; 



(65) 



According to (0), the average is 

(e C(fe z+ (0) (k ) i 1 )C (0) (k , ,t2)) = ^-^^(k + kOe^ 1 * 1 - 421 , (66) 
so integrating ([55] ) over k' gives 

-i J dk(k • B) 2 k c kT x d k C j dt 1 exp {— [— ik • B + \k 2 ](t - ti)| 
•f rft 2 exp{-[-zk-B + AA; 2 ](t-t 2 )}e- 7fcl * 1 ^ 21 . (67) 

In this expression, the integral over t 2 extends from — oo to ti with \t\ — t 2 \ —t\— t 2 , and 
from t\ to t with |ti — t 2 | — t^ — ty- We thus find that 

/ dt 2 exp{-Hk-B + AA; 2 ](t! - t 2 ) -7^ -t 2 |} 



..ifk-B-Afc^ft-ti 1 , 1 



i(k-B-Afc 2 )(i-ii) _| 1 e -7*(*-*i) 



7 | - (-ik • B + AP) 2 -ik • B + A A; 2 - 7fc 



(6* 



so that 

} = Gk-B-Ap ) Gk-B-AP-^ ) ' (69) 
where ( ) refers to the function following dt\ in (|67|). This expression can be written in 

the equivalent form 



dh( )=7 fe ~ 1 



ik-B + Xk 2 ik • B + (Xk 2 + lk ) 



(k • B) 2 + (AP) 2 (k-B) 2 + (Xk 2 + lk f 



(70) 
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so (|67T) becomes 



i / dk(k ■ ~B) 2 k c k~ 1 dkC r f' k 



ik ■ B + Xk 2 ik-B + (Xk 2 + 7 fc ) 



(AA; 2 ) 2 + (k-B) 2 (AA; 2 + 7jt) 2 + (k ■ B) 2 

r 2 



(71) 



We show below that the contributions of the terms of order B in the numerator vanish, so 
(ffl|) becomes 

1 1 



J dk(k ■ Bfk c k- l d k C% 1 



(Xk 2 ) 2 + (k • B) 2 (AP + 7fc ) 2 + (k-B) 2 



Combining (|72]) and fl6"4|) yields the following expression for (|53): 

k-B(AP + 7fc ) 



(v x b) c = -2 J dkk c k~ l d k C 



-|(k-B)V 



(AA; 2 + lk ) 2 + (k • B) 2 
1 



(AA; 2 ) 2 + (k • B) 2 (Xk 2 + lk ) 2 + (k-B) 2 _ 



(72) 



(73) 



To proceed, we adopt a coordinate system (1, 2, 3) such that 



B = (5,0,0) 



has only one component, Bi = B, so that 



k ■ B = k x B . 



Hence 



(vxb) c = 



(AA; 2 ) 2 + (hB) 2 (Xk 2 + lk ) 2 + (k,B) 2 



(74) 



(75) 



(76) 



Note that (v x b) vanishes if -B = 0; as explained earlier, finite B is required to break 
rotational symmetry. 
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In our coordinate system 

dk = k 2 sin 9dpd9dk , (77) 

where 9 = Zj^, and (p is the azimuth in the 2-3 plane. If c = 2 or 3 in (|76|), k c = k sin 9 sin ip 
or k sin 9 cos </?, and as the rest of the integrand is independent of tp, the integral over ip 
vanishes. We are left with c = 1, in which case the integrand is an even function of k\. The 
integral over <p gives 2ir, so 



dk — > —2nk 2 dkd(cos9) = —2nkdk dk\ , 



(7* 



where we have changed variables from (k, 9) to (k, kcos9) = (k, kx). As the integral over k\ 
is from k to — k, and the integrand is even in ki, 



dk — > Airk dk dk\ 



(79) 



with k\ varying from to k. As claimed previously, the terms of order B in (|7T|) would 
have contributed terms of order kf to the integrand in (0); as this is odd in kx, their 
contribution would have vanished when integrated from —k to k. Hence 

r k 1 B(Xk 2 + lk ) 



rco rK 

(vxb), = -2/ Ankdk dk x k x k~ x d k C 

Jk=0 Jki=Q 



(ap + lk ) 2 + (hBy 
1 



(80) 



(Xk 2 ) 2 + {hBf (Xk 2 + 7fc ) 2 + (k l B)\ 
Hence (v x b) is parallel to B as in the result of the classical theory, (|ITD , and the dynamo 
coefficient a is 



a 



r x rk 

-N- / dkd k C / dki 



k=0 



k 2 (Xk 2 + lk ) 



fei=0 



(AA; 2 + lk f + {hBf 
1 



_(Afc 2 ) 2 + (fc 1 S) 2 (Afc 2 + 7fc ) 2 + (fc 1 5) 2 _ 
To carry out the integration over k\ in (|8lD , in the first term in braces we set 

k~B 



AP + 7* 



(82) 



so that 



(A/c 2 + 7fc ) / dkx 



k\\k 2 + lk ) 



/*i=o ± (AP + 7fc ) 2 + (A; lJ B) 2 
In the first term in brackets we set 



.^-tan 1 j 



7] 



kB 
Xk 2 



so that 



_ 2 rk kfdh x x 3 / |7f-7H-tan V 

2 7fc 7o (AP) 2 + (k{B) 2 2lk \ rj 3 

In the second term in brackets, we again use (0), so that 



kfdki 



I^-l t.3 3 



7* * 



/o (A/c 2 + 7fc) 2 + (/ci-B) 2 2 
Combining (§|), (§5]), and (§§), we find that in (gj) 

"7 fc - A/c 2 



k 

fei=0 



where 



so that (PTI) becomes 



a: 



/(0 



7fc + A/c 2 
£ - tan" 1 ^ 



/(O + /fa) 



fc=0 



4irk 2 dk kd k C% 



Ik - Afc 2 
7fe + Afc 2 



/(O + /fa) 



Since the integrand is a function of k alone, 4nk 2 dk = dk and this can be written 

Ik - Afc 2 



a = — dk kd k C'j k 1 



7fc + Afc 2 



/(O + /fa) 



(8 



(8 



(9 



To interpret the quantity kd k C, we multiply the expression on the left hand side of 
), taken for t\ = t, by ik' c to get 



- U. 



>+(o) 



(k)[ik' x z-(°)(k') 



- z 



-(0) 



(k) • w- (0) (k') 



(9 
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where 

w-^(x) eVx z-^(x) . 
From Leslie (1973; eq. 2.16), Q91|) equals 

- h(k)5(k + k') , 

where /i(k) is the transform of the correlation function 



(92) 



(93) 



fc(O = <* +(<9 (x).w-W(x + 



(0), 



From fl5"T| ) with ^ = t, we see that ( |9~T| ) is also equal to 



e cd X i0 \k)ik' c z- {0) (k')) = -ik c {-2ik c k~ l d k C)5{k + k') = -2kd k CS(k + k') . (95) 



(94) 



Equating fl9"3|) with (]95f) , we see that 



so that (p0|) becomes 



| / dkfc(k) 7fc - 1 



^/(0 + /(i) 



(96) 



Ik + A A; 2 

which is our principal result, where h(k),f, £, and 77 are defined by 
(B3) respectively. 



(97) 



(H), (|2J), and 



To compare fl9"T|) with the classical result flUf ), we take the limit — > 0, which 
according to fl82"|) and (|4j) corresponds to £ — ► 0, 77 — > 0, in which case /(£) = f(r)) — > 1/3, 
and the bracket in (|97|) becomes 



3 1 7* + A* 2 , 



2 7fc 
5 7fc + Afc 2 



so 



a (B 0) = -§ / dk/i(k)( 7fc + A A; 



„2\-l 



(9* 



(99) 



-26- 



As the integral of the power spectrum over k is the correlation at zero lag, 
J dkh(k) = (z+W(x, t) ■ V x z-(°)(x, t)) 

; (v (o) +b (o)). Vx (v (o)_ b (o)^ 

r (o) . v x v(°)) - (b(°> • V x b<°>) + (V • (vW x b^)) 

rW-VX v(°l)-(b( '.Vxb<°') (100) 

because ^v^ ) x b^ ^ = according to (2); here we have suppressed the arguments for 
clarity. 

From ( p9|) we conclude that if A is small (but not too small — see below), and in the 
special case that 7^ = 70 is independent of k, 

a(B^0) = -iJdkhCk)^ 1 

V (0) (x,t) • V x v (0) (x,t)) 
|b (0) (x,i)-Vxb (0) (x,t))] . (101) 

The first term agrees with classical result (|T8|) if the damping constant 70 is identified with 
t -1 . The second term has the sign and form expected from the work of Pouquet et al. 
(1976), as indicated in (|T9"|). 
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Our more general result for B — » 0, (99) , differs from the classical result in three ways: 



i) a term proportional to the current helicity is included in the form derived by Pouquet 
et al. (1976); 

ii) allowance is made for variation of the damping constant with k; 

iii) and the effect of finite A is included. 

In the opposite limit B — ► 00, f(u) — ► u~ 2 , so the bracket in ( j97[) becomes 

Ik ~ Xk 2 1 1 lk - Xk 2 fa + Afc 2 ) 2 (AA: 2 ) 2 _ ( lk\ 2 ( , 

~/ k + \k 2 u 2 v 2 j k + Xk 2 (kBf (kBf \kBJ ' 1 ] 
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independent of A. Thus, 



a(B — > oo) 



2B 



dkh(k)k- 2 lk . 



(103) 



An inverse dependence on B was cited by Krause and Radler (1980), but with a different 
factor. 

Finally, we consider the simplification introduced when the spectral density C(k) is 
effectively concentrated at some wave number h*, which if k* ~ 2k is a crude representation 



of Figure 1 of the numerical result of Pouquet et al. (1976), where E\ peaks at A; ~ 1, E t 
peaks at k = 3, and the energy beyond k = 3 is Alfvenic, so it does not contribute to a. 
Then from (H), @), and ©, 

7* - A/c 2 



where 7* = 7fc(&*); 



and 



7* + A A;; 

KB 
7* + AA; 2 

kSB 



-/(£*) + f(v*) 



(104) 



(105) 



(106) 



In the light of (|101|), then, 

"V ).Vxv(°)\-(b^-Vxb( ' 



a 



2 '* 



7* - AA;: 



/(£*) + 



(107) 



7* + A A; 2 

so the dependence upon B is entirely in the multiplicative factor. If we define a magnetic 
Reynolds number by 



R 



M 



7* 
\kl 



(108) 



and a nondimensional mean magnetic field by 



B 

(7* A*) 



(109) 
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we can write the bracket in ( 107 ) in the form, 

F(R M , P) = f^/ (P^rr) + fiPRu) ■ (no) 

As /3 — > 0, this approaches 



F(R M ,(3^ 0) = | + l) = §—£-73 , (HI) 

6 \ R M + 1 / 7* + Xki 

which, when inserted into (|107|) , yields agreement with ( |101|) , our previous result for B — > 0. 

However, one must be careful about this procedure if Rm is ^> 1, as is usually the 
case, for even if (3 <C 1, (3 may be ^> R M , so that the second term on the right of ( |110| ) is 
effectively rather than 1/3. In this case 

F(R M ^oo,P^0) = l, (112) 

and a is given by 1/2 the value in ( |101|) . Evidently, the classical expression for a (as 
modified by Pouquet et al.) is off by a factor of 2 if Rm ^> 1. Figure 1 is a plot of 
F(R M ,/3), compared with F*(R M ,f3) = (1 + Rm/3 2 )' 1 from Cattaneo and Hughes (1996). 
As F*(R M ,P) was not proposed to be accurate for R M < 1, the curves of F* for R M < 1 
should be ignored. In qualitative agreement with F(Rm,P), F* decreases as Rm increases, 
but of critical importance, it decreases to 0, rather than remaining finite as Rm — * oo, the 
case for the galactic dynamo. 



4. Discussion 

Our general result for a, (|9T|), displays significant features of the a-effect due to driven 
helical turbulence. Our discussion of ( |101| ), above, makes clear that the values of v and b 
referred to in h(k) are v 1 - ) and h^°\ the values of the small-scale velocity and magnetic 
fields that apply if B = 0. 
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As indicated in the text, it is in principle possible to calculate v^ ) and tr ', and hence 
h(k), from numerical simulations of driven helical turbulence with B = 0. Then evaluating 
a from ( |97|) requires the value of 7^, which can be calculated from ( |6lD in terms of the 
statistics of the zero-order state. One should also calculate a directly from its definition 
using simulations with various values of B, and then compare the results with (|97|). 

Another feature of our work is that the contribution proportional to the current 
helicity, • V x b^ ^, first described by Pouquet et al. (1976), emerges naturally from 
our work. To give it a physical interpretation, we refer to a derivation we carried out in 
terms of the variables v^ ) and b^ ) to check the derivation presented here in terms of z^°\ 
Referring to (|25J), we find that at first order in B, the current helicity term is traceable 
to X b(°)^}, as opposed to (v^ x b^ 1 ^, which gives rise to the kinetic helicity term, 
where the velocity perturbation is due to the first-order Lorentz force (V x b^ ^) x B 
as explained by Pouquet et al. (1976; p. 332). 

Current helicity is discussed in the Soviet and Russian literature. Vainshtein (1972) and 
Vainshtein and Zeldovich (1972) called attention to effects of current helicity (which they 
referred to as "magnetic gyrotropy"). They argued on physical grounds that as B grows, it 
induces current helicity, which, because its sign is opposite to that of (v ■ V x v), reduces a, 
causing dynamo activity to cease at some finite value of B. As we have explained, Pouquet 
et al. (1976) showed that current helicity is important at all values of B, and on the basis 
of their spectral equations, found that the effect described by Vainshtein (1972) would be 
smaller and of the opposite sign than he claimed. Vainshtein and Kichanitov (1983) accept 
the interpretation of Pouquet et al. (1976). They also derive a differential equation for the 
small-scale current helicity from the conservation of magnetic helicity as did Kleeorin and 
Ruzmaikin et al. (1982); see Kleeorin et al. (1995). We discuss these matters further in a 
forthcoming paper (Blackman and Field 1998). 
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The present work relies upon r <C t ed to evaluate the perturbations of v and b due to 
B (although not in evaluating v^ ^ and b( )). This is a reasonable approximation because 
the dominant contribution to a comes from wavenumbers < 3k , which are relatively 
unaffected by magnetic back reaction, and therefore should follow Pope's (1994) finding 
that in hydrodynamic turbulence, r < t e d- Unfortunately, the inequality is not as strong as 
we desire. 



5. Conclusions 

We have found an analytic formula for the dynamo coefficient a based upon a 
perturbation expansion in B, the magnitude of the large-scale magnetic field, which is 
responsible for perturbing isotropic helical turbulence in such a way as to produce a 
turbulent emf along B. Our formula gives a in terms of B and the magnetic diffusivity 
A, together with the spectra of kinetic and current helicities, and a damping coefficient. 
Both of the latter are calculable from a simulation of incompressible isotropic helical MHD 
turbulence. For small values of B, our results agree with the classical results as modified 
by Pouquet et al. (1976). For large B, we find that two terms contribute to a. The first 
term is independent of the magnetic diffusivity A, and is similar to the expression proposed 
by Kraichnan (1979). The second term, which depends on A in a way reminiscent of the 
expressions suggested by Cattaneo and his collaborators, vanishes in the limit A — > 0, as 
predicted by Cattaneo, et al, but contrary to Cattaneo et al, the first term remains finite, 
and so a is reduced only by a factor of 2 for large values of B. 

We are grateful to Amitava Bhattacharjee, B. D. G. Chandran, and Eugene Parker for 
their insightful comments, and to the Aspen Center for Physics, where part of the work was 
done. We are also grateful to Mr. Dimitri Kalaitzidakis for his hospitality and allowing us 
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A. Energy Cascade 

Start with the Navier-Stokes equation, which describes incompressible hydrodynamic 
(HD) turbulence: 

a (V = -vVv-Vp + i/V 2 v + f, (Al) 
where p = pressure -hp. Define the kinetic energy E v by 

E v = (\v 2 ) . (A2) 

Then the scalar product of (Al) with v shows that 

d t E v = -u(u 2 ) + (f -v) , (A3) 

where 

(jj = V x v (A4) 

is the vorticity. Hence in unforced ideal HD, E v is conserved. The only effect of the 
nonlinear term — v ■ Vv in (A.l) is to redistribute E^ , the kinetic energy spectrum, over k, 
keeping E v = j E\dk constant. Numerical simulations show that this redistribution takes 
the form of a direct cascade, that is, from lower k to higher k. 

Now add a magnetic field B, so the momentum equation becomes (fZ5|): 

d t v = -v ■ Vv - Vp + B • VB - V\B 2 + z/V 2 v + f 

= -v- Vv + B- VB- VP + z/V 2 v + f , (A5) 

where 

P = p+ i B 2 . (A6) 
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The equation for E v becomes 

d t E v = (v • (B • V)B) - v L 2 ) + (f • v) , (A7) 



so E v is not conserved in unforced ideal MHD because the first term allows magnetic 
energy to exchange with kinetic energy. To account for this, we write the induction equation 
governing B: 

d t B = -v • VB + B • Vv + z/V 2 B , (A8) 

where we have specialized to the case of unit magnetic Prandtl number, \jv. From the 
scalar product of (A. 8) with B we obtain an equation for the magnetic energy E M = (i£ 2 ), 

d t E M = (B (B V)v) - v ( J 2 ) , (A9) 

where 

J = V x B . (A10) 

E M is not conserved in unforced ideal MHD because of exchange with E v . However, we 
note that in (A. 10) 

(B • (B • V)v) = (B • V(v • B)) - (v • (B • V)B) 
= (V • [B(v • B)]) - (v • (B • V)B) 
= -(v(B-V)B) (All) 

since spatial averaging converts the divergence to a vanishing surface integral. Hence 

d t E M = - (v • (B • V)B) - v (j 2 ) . (A12) 



The energy-exchange term in (A12) is the negative of that in (A7), indicating that energy 
gained by E v is lost by E M and vice-versa. Hence if v — f = 0, the sum of (A. 7) and (A. 13) 
yields 

d t E = d t (E v + E M ) = (A13) 
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for the conservation of the total energy E in unforced ideal MHD. By analogy with E in 
HD, the second nonlinear terms in (A. 5) and (A. 8) for MHD can only redistribute Ef. over 
k. Numerical studies like those of Pouquet et al. (1976) show that in the non helical case, 
the effect of the nonlinear terms in MHD is a direct cascade of total energy. However, 
because E v and E M are not separately conserved, it is not obvious that both E v and E M 
individually cascade directly in general. 

To answer this question, we turn to the Elssaser variables for the small-scale field z^, 
which in the case A = v satisfy fl35"[) : 



d t z ± = -z 



T • Vz± - 



VP + z/V 2 z ± + f ± B • Vz 



± 



(A14) 



Define 




(A15) 



Then the scalar product of (A. 15) with z 1 * 1 yields 







(A16) 



Since V • z 



± 



0, this gives 



d t E 



(V-^z^) 2 ])-^-^)) 

v (z ± ■ V x w ±) + (f • z ± ) ± (V ■ [Bi(^) 2 ] ) 



(A17) 



where 



c; ± = V x z 



± 



w± j • 



(A18) 
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Since 

z ± • V x w ± = V'f^xzj+Z-Vxz 1 

= V-^xz^ + ^f (A19) 

and averages of divergences vanish, 

^ ± = -z/((cu ± )) 2 + (f •z ± ) , (A20) 

an equation identical in form to (A. 3). We conclude that if A = v both E + and E~ are 
individually conserved in unforced ideal MHD, and that the effect of nonlinear interactions 
is to redistribute the two energy spectra E£ and ££" individually. Note that the presence 
or absence of a mean magnetic field B does not affect this conclusion. 

That leaves the question as to whether the redistribution is a direct cascade. Since 

E ± = (|(, ± ) 2 ) = (|(v±B) 2 ) 

= (y) + (^ 2 )±(v-B) 

= E v + E M ± 2K = E ± 2K , (A21) 

where 

K = ±(v.B) (A22) 

is the cross helicity (Biskamp 1993, p. 179), it follows from the conservation of E + , E~ , and 
E that K is also conserved in unforced ideal MHD if A = v. As pointed out by Pouquet 
et al. (1976), if the initial state is assumed to be statistically invariant under b — > — b, it 
vanishes initially and thus remains equal to zero, and from (A21) we have 

E + = E~ = E . (A23) 

As E cascades directly, it follows that E + and E~ cascade directly if A = v. 



B. Correlations 

Here we show that (|58|) is a reasonable representation. Let 

V=(z + (°)(k,t)xz-(°)(k',t 1 )) , (Bl) 
so that the average in fl57|) is V c . From ( |50D 



V = (2vr)- 6 / / cfo«/xV i(k - x + k ' x) (z + W(x,0) x z-(°)(x',t 



(2tt)- 6 / / rfxdre- i(k+k ' )x - tk '- r (z + (°)(0,0) x Z -W(r,r)) , (B2) 



where we have put t\ = t + r and x' = x + r and have used stationarity and homogeneity. 
Carrying out the integration over x gives 

V = (27r)- 3 <5(k + k')i(k', r) , (B3) 



where 



If we write 



(k', t)= J dre- ik ' r (z +(0) (0, 0) x Z - (0) (r, r)) . (B4) 

Q(r,r) = (v(0,0)x v(r,r)> , (B5) 

S(r,r) = (b(0,0)xb(r,r)> , (B6) 

P„(r,r) = (v(0,0)xb(r,r)> , (B7) 



and 



P 6 (v,r) = (b(0,0)xv(r,r)) , (B8) 

where we have omitted superscripts for clarity, then the average in (B4) is, from (p0|) and 
(0), 

z+(°)(0,0) x z-(°)(r,r)) = Q(r, r) - S(r, r) 

-P„(r,r) + P 6 (r,r) , (B9) 



so 

i = Q(k', r) - S(k', r) - P w (k', r) + P 6 (k', r) , (BIO) 

where 

Q(k',r) = |rfre- ik ' v Q(r,r), (Bll) 

etc. One of us has shown (Chou & Fish 1998) that the terms P„ and P 6 do not contribute 
to the final result, so we ignore them in what follows. 

Consider 

Q(r,r) = (v(0,0)xv(r,r)> 

= (v(-r,-r)xv(0,0)> 

= _( v (0,0)xv(-r,-r)) , (B12) 

where we have used homogeneity and stationarity and flipped the cross product. Thus 

Q(k', T ) = - J dre- ik ' r (v(0, 0) x v(-r, -r)) 
= - J rfre ik '- p (v(0, 0) x v(p, -r)) 

= -Q(-k',-r). (B13) 

Krause and Radler (1980, p. 75) show that quantities like Q are odd functions of k', so 
(B13) implies that 

Q(k',r) = Q(k',-r). (B14) 
Since S(r, r) has a similar structure, 

S(k',r) = S(k',-r), (B15) 

so that I(k',r) is an even function of r, or, equivalently, a function of |r|, which we shall 
denote (27r) 3 R(k', |r|). Hence from (B3) 

V c = 5(k + k')R c (k',\r\) . (B16) 
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Because R is composed of cross products, we introduce a tensor R de (k', |r|) such that 

R c (k',\r\) = e^it^k', |r|) , (B17) 

so that 

V c = e cde 5(k + k')R de (k\\r\) . (B18) 

The last step is to assume that Rde(k', \t\) is a function of k' times a function of |r| 
(which may depend upon k — |k| = |k'|): 

R de (k',\r\) = R de (k')f k (\r\) , (B19) 

where by taking /fe(0) = 1, we assure that R de (k') represents the maximum value of the 
correlation. As a natural choice for /fc(|r|) we take 

/*(M) = e-^W , (B20) 

where 7^ is a fc-dependent inverse correlation time, in the spirit of the eddy-damping 
approximation of Pouquet et al. (1976) and Kraichnan (1979). Then 

V c = e cde 5(k + k , ) J R ( i e (k')e-^l T l , (B21) 

which is fl58|). 
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Fig. 1. — Three functions of the magnetic Reynolds number R M and the normalized magnetic 
field f3 = B/vo are plotted and compared. F(Rm,/3), defined in (110) and represented by 
solid lines, is our result for the dependence of a on R M for (3 = 10~ 5 , 3 x 1CT 3 , and 1. 
F(R M) f3 — > 0), defined in (109) and represented by the circles, is our result for f3 — > 0. 
F*(R M) f3) = (1 + i?M/9 2 ) -1 is the dependence of a suggested by Cattaneo & Hughes (1996); 
as it was not intended to be accurate for R M < 1, points for R M < 1 should be ignored. 



